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ABSTRACT 





A mathematical model of a branching stochastic process utilizing 


| 


generating functions is presented, The probability distribution 
of the number of members of the process at discrete time periods, 


2, and the probability of extinction is discussed, When there is 


n? 
a non-zero probability of surviving indefinitely, the normed random 
variables, 2,/E(zn), converge with probability one; the cumulative 
distribution of this random variable is absolutely continuous. The 
time until extinction and the total number of members of the process 
is examined when the probability of extinction is one, The distri- 
bution of the 2, given that z, is not zero is discussed for this case, 
The maximum likelihood estimates for the probabilities involved in 
the process are determined, An example is given of a branching pro- 


cess in which the probabilities are dependent on time and a solution 


is found using Laplace transform methods, 
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NOMENCLATURE 


Meaning 

the probability of the event A 

the conditional probability of A given B 
the joint probability of A and B 

the number of members of the rth generation 


th member of the 


the number of offspring of the i 
nth generation 

the expected value of 2) 

the variance of 24 

generating function of Zy 
probability of r members in the first generation 
generating function of 2, 

probability of r members in the n&) generation 
probability of extinction 

a normed random variable (“Ay 

the limit of w, as n approaches infinity 

the cumulative probability function of w,, 

the characteristic function of w,, 

asymptotic distribution branching from f(s) 
probability generating function 

the asymptctic distrieetion branching from k(s) 
the characteristic function of H(u) 

probability that 2, 4 0 


probability that 2, equals r given that it is not 


equal to zero 











Yn 


F,(s) 


number of members of the n&h generation given that it is 
not zero 

probability generating function of y,, 

number of generations to extinction 

moment generating function of N 

total number of individuals in the process 

probability generating function of S 

number of individuals in the m*) generation who have 
exactly r descendents in the (m 1) st generation 


total number of individuals in the first n generations 


vi 











INTRODUCTION 


The present paper is an essay on the theory and application of 
discrete branching stochastic processes, Most of the esymptotic 
theory was extracted from the work of Harris [10] who derived many 
of the relationships for this class of stochastic processes, Let 
us introduce these processes by means of a historical example, 

In 187), one of the first investigations into discrete branching 
processes was made by the Rev. H. W. Watson. Dr. F, Galton had 
posed the problem of extinction of family names to the mathematical 
community in England since many names, prominent in history, were 
no longer to be found. The problem was posed as follows: Assume 
a large nation with N adult males, each with a differe:t surname, 
What proportion of names would be extinct after r generations? In 
the rth generation, how many people would have the same name? Rev. 
watson [8] derived some relationships that enabled him to compute 
the probabality of extinction and made some sample calculations to 
demonstrate the difficulty of computing numerical answers, Notable 
among these relatioriships was the functional relation of the proba- 
bility generating functions involved (section III, theorem 1). 

Branching processes received some attention thereafter ard the 
theory was applied to problems of gene mutation [7], population 
growth [12], cosmic radiation [1], ionization fluctuations [18], 
neutron production in atomic piles bil and cascade theory [20]. 

In this paper we will review the theory of simple discrete 
branching processes and some of its applications, Let us consider 
the original problem of survival of family names ard state the 
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basic problem of branching processes in that framework, Suppose 
an immigrant by the name of Door comes to the United States and 
there is rio one else here having that surname, What is the prob- 
ability that the name survives r generations and how many men will 
there be having that last name? 

There is a probability that Mr, Door will either not marry or 
not have any male children; there are other probabilities associated 
with the number of sons he may have, We will aot these probabili- 
ties by py, 1.e., the probability of no sons will be po, of r sons, 
Pre We need a symbol for the number of individuals with the name 
of Door in the nt) generation, so let 2,, be the number of individuals 
with the name of Door in the nth gereration, Then we can write 
P(z) = 1) = py indicating that the probability of r individuals with 
the name of Door in the first generation is py. If Mr. Door coes 
have sons, each of them may have children and tre process is continued. 
This may be represented as a tree type graph with an initial node 
representing his sons, and so on, A method of analysis of the 
problem based on graph theory was proposed by Otter [19], In this 
paper the method of generating functions is used, 

We will assume that the number of sons born to one of Mr. 

Door's sons follows the same probability law as that of Mr. Door. 
In other words, the probability that Door, jr., has n children is 
Bye linen, if Ke) is the number of offspring of the 4th individual 


m 
in the nth generation, we have P(zn+] = r|Zn = m) = P(= x{n) =n). 


This is the outline of the problem, the notation, and the assump- 
tions made in this paper. Before taking up the mathematical theory, 


some other problems will be posed in this framework. 
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In atomic piles, the number of free neutrons is an important 
factor, If we start with one neutron and a fission takes place, the 
probability that r neutrons are released is a nuclear constant which 
must be empirically determined for the amount and arrangement of active 
material, The probability of no fission (no descendents) is the 
average probability of leakage or absorption. Am atomic pile in 
a nuclear reactor is subcritical when the reactor is not operating 
and there is no possibility of an explosion, emithe probability 
that the number of fissions increases sufficiently to cause an 
explosion is zero; the mean of the probability law for a subcritical 
mass is less than one. As the amount of fissionable material is 
increased, the probability of an explosion finally becomes greater 
than zero. When this occurs, the mass is termed supercritical and 
the mean of the Pecancvestiean law for this situation is soumacl to be 
greater than one, A nuclear reactor is termed critical when it is , 
operating; the mean of its probability distribution is exactly one. 

The organization of the paper is as follows: The general math- 
ematical character of a discrete branching process is presented in 
section I. It will be seen that the mean number of progeny per 
eet plays an important role in the asymptotic theory, This 
theory is presented in section II for the mean of the probability 
distribution greater than one and in section III if it is ‘not 
greater than one, The maximum likelihood estimation of the distri- 
bution of progeny per individual is given in section IV. An example 
of an ionization cascade is presented in -sectiomw Yow Its asta iene 
general branching process than those treated in the model used in the 


earlier sections and its treatment suggests now the methods of Markov 


processes may be applied. 
<o= 








GENERAL 
The structure of branching processes was presented in the 
introduction where several assumptions were made ard some notation 
was introduced, It is appropriate tow to summarize the assumptions 
made in this paper and to extend the notation, 
Consider two related classes of random variables: 
Z, ‘the number of "individuals" in the nth generation 
xin) the number of "sons" of the i? individual in the 
nth generation, 
They are related by the equations: 
Zine! (n-I) 


Zhe G Pe LOL =e cy ges fo te 


The quantity 2,5 is the number of the initial group who start the 
process; Z) is assumed to be one, since this does not affect the 
generality of the theory as is pointed out at the end of this section. 
The probability distribution of the first generation is given by p,, 
1.€., the probability that there are r members of the first generation 
1S p,- 

The following assumptions are made throughout this paper: 
(1) Given z,, is equal to m, the random variable, ans] , is the sum 


of m independent random variables, each one of which has the same 


distribution as 2) : Py.” = r) =r 


(2) The variance of 2) Laetinive., 


(3) p,# 1, for all r ; otherwise the process degenerates to a simple 


arithmetic calculation. 
(4) po+ py 41, otherwise we have the binomial case whose treatment 
“in the context of branching processes is trivial, 


oie 








The following notation is used throughout: 


oO 


JP =El2}=S er 


the mean of z 
r=0 } 1 


Za 
6* = Var (2,) = ir Pe) the variance of 2) 


Cc 
= r 
F(s) = = P ry the probability generating function 


of 2, (s is a complex variable) 


car” mics. =/) 


- | 

Fn (S) ie oe the probability generating function 
of 2. 

From the assumptions listed above, it is clear that : 
nats) = 8; a ee f'(s) and f"(s) are continuous on the set of 
points consisting of |s|<1 and s=1. 

Several authors (see [2], [9], [21]) have examined the generating 
functions and have found some basic relationships. Some of these are 
presented below, 


iMieerem 2. £,,(s)-is the nth functional iterate of f(s): 
FRi(s}=F (Fr(5)) = Frey (F (SI) 


Proof: From the initial assumptions, the re are independent for all 


(n) Zn-1 (nel) 
teane 1, P(X; = K) =P, ZnS x; 


It follows that Cierra ict is the coefficient of sk 


in the power series expansion of (f(s))2 since the generating function 
of the sum of j mutually independent random variables is the product 
of their generating functions. See Feller [6, p. 250ff.]. 

ie 








Hence, F, (S) = - Pj (c (s))” = F (€(s)} 


f(s) “ZB (eis) = FaCFts)) = FCE(F(S32 = F Le tst) 


i 60 
Frisk= = py Cris = Fpl FSi) = Flt. (FLESH)... .1) 


The coefficient of s® in tcs) is the probability that there 
are r individuals in the nth generation, Thus the questions posed 
by Dr. Galton may be answered in the framework of the present model. 
The constant term, p,,, is) the probability that the family name is 
extinct by the nt gereration, | 


The mean and the variance of Z,, car be determined using theorem l. 


E (Zn) =n" Var (7n) = SMe yealh fé¥} 


Theorem 2, 


~ nom f>s 
Proofs (i) = > rp = E (2h) 
(i) = = ee p= E (Zn) = E (Zp} 


Therefore, Vor (Zn) = Fi) i Ft) = C6 "(]”, 


By theorem 1 mre CS r EE USa es ECF si] 


Taking derivatives, Se (s) = F, CF (si f(s) = F CF, (s)] , (s! 


= (Soe 





E (s) = F" C€(s)] Peter + E Tesi e “(s) 


FT EFR(SIICEIS™ + €'CEA(SIT (5) 
Hence heat lve! E (Zp) sy 


F, CRMC EOI +f, Cee) = EULER WILE! (o> + CE, (IIE C0 





Solving for f!! (1) gives 


Set ygy = Ear WW 6.0) EW 
i ade - Foil) 


_ (or =u + U2) | Te A) 


Jam = fe 
a,,n n 
= 6 ML (M -|) _ n an 
Mey Via T/A 
Thus Var(z,) = on 0 (0-1) or 
ius Var(z,,) jada pa i jAFi 


For the case Mal , L'Hopital's Rule can be used: 


ee Zn an-| = ra nel) 
aM -| 


i 6*(Qn-n) 


a 


a =| 





jae 


Sicmmavat CZ») — ne for Mel 


Higher moments, if they exist, may be found by a similar pro- 
cess utilizing the higher derivatives of fe (s)- 

A subject of great interest in any branching process is the 
probability of extinction of the species.: The Fundamental Theorem 


re 








of branching processes tells us that this probability is a fixed 
point of the generating function. The proof given follows that of 
Bharucha-Reid [2]. First, we note that if "a" is the probability of 


ee “14 - 13 
xtinction, then a Lim Das Lim, EI ES)) 


Theorem 3. The probability of extinction, a, is the smallest pos- 


iawemnurnber, p, such that : 


p=Fipl= > BP 


Further, a = 1 if and only if jf =], 
Proof: We see that a = 0 if and only if py = 0 and f(1) = SO 
that Ofae 4 as it should be. 


“To show that a = f(a), we note that since p,,, = f,(0), 


= FCF,(003 


Prato 


SO a = Lim Pa Soe ti (0) = lim tie aC) “ie £(Pn.1,0) = tla) 


Now f(s) is a power series in s with positive coefficients and hence 
strictly increasing for real arguments, If O<a<b, then f(a) <f(b). 
Also, since py,20 for all n, then azZ0O. 


Let p20 be such that p= f(p). It follows that 


B= Flo) <Fl(p)=p 
Assuming p,o€ p, then f,(0)< p and 


nts = Fae (0) =F LF niol) € Fp) = p 


So, by induction, pp o¢p for all n. 
Then lim €p and afp. 
It remains to show that a = 1 if and only if 703 ae 


ie 





Assume first that a = 1 and show that /4 <j 

Now, it may be noted that: 
f' (s) 2 0, f"(s)Z0, so f(s) is convex in (0,1) 
Hi iGo) we? Fo} 10) ws) aaa 


f(s)#s for s in (0,1) since 1 isthe minimum positive number 


satisfying s = f(s). ELS Ry 
AY 
Therefore, 1 - f(s)<1-+s and. { 
a) is bounded by 1 and is 
monotonically increasing e 
with s for s in (0,1). 
i S 


Therefore, f'(1) exists, f'(1)é1 and Med, 

Now assume M€4, Then, since f(s) is convex, f'(s) is either 
constant in (O51) or Serictie increasing with s, By assumption, 
p, #1 sop >0 and f'(s)< 1 in either case for s in (0,1). 


i 
By the mean value theorem, F\\l- F(s)= F (x){i-s) where s<xel 


» 


Then, |-F(s) ='(x)(I-s] and 
I- Fis) < [-s 
Fis}>S - for ofsél 


Wienewore, a= i, 

foe eoren 3, it is a simple matter to determine the limits 
Of bhegp eetior 1x 0x 
Theorem . A discrete pronenine process either dies out or becomes 


extremely large, i.e., Jim pers Out omuiets0e 


aon 





Proof: Let t be a real number such that 0< t<1 and recall the 
iteration property of f,(t). 

Case 1.44 €] Then f(t) 7t (see proof of theorem 3), and 
mgm f5(t)e. . oni f,(b)< . « « El. 
Therefore, tim f(t) = 1s a, 


Case II.uzvl and t<a (p,70). Then f(t)>t since f(s) is 





convex and f(t)Z f,(t)< oe o SL (t)K o o @ aa Ao 





em 


Therefore, iim, f(t) =e 


Case IlI.u>i and ta. 
Then f(t)Zt due to the convexity of f(s). 


femeer f(t) >f5(t)r~... pi,(t)7 oe o> a. Pp 





Therefore, tim f(t) = a, 4. S 
oo oo 
Now f(t) = 2 Pip t” Spent 2. p,-t by definition, Taking 


the limit of the last expression gives: 

of 
lager (G) = lim lim tT =a from the above three cases, 
n»o n't) JB Poot 133 2 Prr 


. a ers) | 
Theref Is r =0, Since the series converges uniforml 
Ret One auth Part e) in seri O g y 


7h) 


in an interval one may pass to the limit before summing, 
1°, @) 
Therefore 2 lim Ppt? © 0 and lim ppp = 0 for r f 0 by the 
uniqueness of power series, 
Theorens Beand J give us some interesting information. If we 
are working with an atomic pile, the probability is one that it will 


die out unless it is suvercriticals; then the probability that it 


will blow up is l = a, 


so 








Lotka [17] found that the probability law of the number of sons 
aman might have was closely approximated by a type of geometric 
distribution, namely, p, = 0.4825, Pp, = (0,212€) (0,5893)K-1 kZ1 


based on a 1920 census, ‘ It follows that: 


Reel yaaa 7 
f(s) = 0.4625 + OsereS fae P's) a 1. 2622 
Solving f(s) = s, it is found that a = 0.6226, Based on these values 
the probability of extinction of a family name is (0.6226)"™ where 
nis the number of male members of the family at the present time. 
This formula follows from the fact thar if there are r members of 
the zeroth generation then the generating function for the first 
pemeration is (f(s))*, Then f,(s) = f,.,((f(s))*); af f(a) =a, then 


(f(a))? = a? and a¥ is the probability of extinction. 


=lle 





MEAN GREATER THAN ONE 


In this section we will consider the case where is strictly 
greater than one and the family has a non-zero probability of 
surviving indefinitely, Since we have determined that z,, may 
become extremely large, we will work with a normed random variable, 


Z | 
n 
Let W, = —w- . The mean and variance of w, are easily deter- 





mined, 
Zn x . 
Elwn)= E (Zor or E (zp) = 1 
_% a, N{,,n 
ull Zia eel gg Poy an 
me (Wr } = E (7a = she ~~ + pA | 
ol | 
= za55—— | ~ jan + 1 
4 a ape 
Var (Wp] ope) (| wi 


Treatment of the convergence of w, may be found in Harris | 10 | 
and Bharucha-Reid [2]. We will prove first, the convergence in 
probability and then the convergence with probability one. 

Theorem 5. If u>i, the random variables, w,, converge in probability 
to a random variable w. 


Proof, Let n, m, be integers and n>Mm, 


Me 


> | [coefficient of sd in(fy__(s))™] 


-_= 
~ 


| oa 


o_ 


( From (S))° 


$s 


| 
Q_ 
a) 


ee 





i} 


E | ele een) r[ Fem (sy FremlS |. 


rol me 
r (1) ail 


RE 


~ 


nue 





EC (Wy-wml@] = E (we)- E(w) 








lm El Wy -Win)” = 0 
Mm, Ne 
Hence the random variables converge in mean Square and by a theorem 
of Kolmogorov [13, p.3l, I] , in probability. 


Convergence with probability oné may be proven with the aid of 


Doob's [|] Martingale convergence theorem. 


Sige 





A discrete parameter martingale is a stochastic process such that; 

a) E (ly ni) ¢ od for all n 

b) EB (yyy | Yns Yn-ls © © © 5 ¥1) ® In 
The theorem may be stated as follows, Let (y,, n2zQ) be a discrete 
parameter martingale, Then E(\y,|)€ E(\yi|)$E( yo )Z. © » 
If lim aa =k Zo, then iim Yn = y” exixts with probability one 
and E(\y*|) < 


Considering our random variable wy, , E(w,) = 1 for all n, 


nth 7 

= | | = a 

E Wha | We) = E | sarin ad 
7a 7 


fn 


pr 
W, 


if 


Hence the (w,) form a Ciscret* parameter martingale and we have 
Theorem ha, If uw>1, the random variables, w, , converge with probability 
one to a random variable w. 

Since w, converges to w, we may discuss the probability distri- 


bution of w. 
; a) 
Let G,(u) = P(wyéu), — * E( e“n”) a jes dc(u), 
G(u) = P(wu), P(s) = ie eSU dG(u). The distribution G(u) can be 
0 

called the asymptotic distribution branching from f(s). There is a 
very interesting functional relation between the moment generating 
function (0(s) of w and the probability generating function f(s) as 


given by the following theorem, 


ie 





theoren 5. (P(us) = F (UIs) 


Spoor nls) 


Hf 


={,4hnS)_¢{, 25) << 
Ele J= Ele) => ep expla $) 


if 


F,(e x) 


Oris) = Eny, (em) = ELE (eH) 


7 FE, (s)} 


lim Cn (s) = Gls ,therefore, 
n-»cd ‘ 
SU 
lin Yr, | vA r L dG, (v 


aiid Giv) = (p(s) 


Upon taking the limit of Pry (Ms)= FC Chaa (s)) the desired result 
is obtained, 

It also follows that the k*® moment of w exists and is finite 
if and only if the kt® moment of 2) exists and is finite, 

Since G n(0) = P(Wwn¢ 0) = py , then Lim, G,(0) = a and P(w = O)z2a, 
The following theorem will be used later, 
Theorem 6. Let G,(u), Go(u) be distributions having equal first 
moments and finite second momerits such that their characteristic 


functions 0, (it) and (Po (it) satisfy Lp, (itu)= Fl LY, [core = 1,2 


Then G 1) = Go(u) 


ice 





K-I | : 
proof: (p(it) = +2 = (it) o(t") 


= 


for |t| small by a special form of McLaurin's theorem for small 


values of |t]. See[Cramer 3, p.27] . Then, 


| os 
[Ulitl- Galiebe |} 4 4 Se ee ey" + O1L ee 
Ns! 
a (2) 
A a (us it)" + 0 it where ip (st) se ed 
=) ) 


1, (Lt) = We (et)| 


) 
(oe oats SSAA. ac 5 19 


te js (or ov) + 0 (2) 


i 


teretore, | Y,(6t)- Ua (it)l= t” alt) where lim @lt) <0 


Then | Up, (LEAL) ~ Lp. (ut) | = | FY, (ut = FOU, (.#))| by theorem 5 


& MANY lu, )-U, | (re)! since [f (s)]éu when {5/2 | 
Hence eB (it) ae 
and pa] G(mmt)| < 16 /¢)| implying that 14 1@(0)| € | Alo! 


But pA >1 so 4'(t) must be identically zero, 


Therefore | (ut) = Yet)! =O 
and UD, | (ut) = = Ups (ct) 


Finally G,(u) = Go(u) by the uniqueness of moment generating functions. 
In order to determine some prcperties of G(u), it will be con- 
venient to work with a closely related function H(u) which is defined 


as the asymptotic distribution branching from k(s) where 
kis) = £osli-altal -G_ 
I-O 
allée 








To demonstrate that k(s) is a generating function use the 


fact that sav] implies f(a) = a¢1. Then 


. Flaj-~ a u 
OS ag, = 0 
Ie {4} = Hine = | 
| -O. 


Clearly k(s) is monotonically increasing since, if ts, then 


Ct li-a)taj-a ((slt-altal— 
Pi) ede toltaj~o va [stt-alt aj ~q 2) 


l-a a {=Ch 


Therefore k(s) is a probability generating function, 


It is interesting to note the following properties. 


Kls\< €[sll-a}+al (1-a) c 


c' (i-aQlhtcG 
C'CS i 


£2 oe 
= _! SS y Se r 
Kis) = --e-( > pist-aital ; wae 
r=0 Ys 
iWkEn) 
ies k \Ked 7a Bee ae avala far ee 
The coefficient of sk is (1-a)“"(p +p, Kar. ..th wy iP Got..e! 
and 4=PRtepne - 3 Pose > 3 = f'(a), 


a 
t ! = vs > B ile f ar it Similar t 
Let Yls) = S. 2 of bt vu) y a line o gument similar to 
the proof of theorem 5, it can be shown that W (eS) =KL Hs 


ive: 


/ 





To demonstrate that k{s) is a generating function use the 


fact that a>] implies f(a) = a¢1. Then 





flal- GQ | 
K (o)= a =O 


efi) > tia 
1-A 


fi 
_e 


Clearly k(s) is monotonically increasing since, 


Kit}e= fitu-oltaj-a | f(stt- altal—q A Kis] 
ha a =a 


if tees, ven 


Therefore k(s) is a probability generating function. 


It is interesting to note the following properties. 


K'(s)= Bielnclzai(\-a) = f'(slt-alta] 
Kir) = ¢' ty) = fh 
Ke is)= €"Csti-aj+a](l-al 
K' (j} = oa I< G] 
KU) EK O)- CK ONT = ot (l-a)-a leu) 
09. 8 

= } _ | at re iv 

K{s] Tal 2. &iSikaj+a} zi a,§ 


k-1 KEN) r 
The coefficient of s* is (l-a) (Pp + ron KA +. + Oo $ ee 


K i fcen i 
and qG,= P, tlpa + 3Rat en Smee {e), 


Let W/s$) = Se Sd) bt wu) 


the proof of theorem 5, it can be shown that YW Ue 5) =KL[ Ws)). 


» By a line of argument similar to 


Sige 





We now proceed to determine a representation for the random 


variable w. First we determine the relationship between |) (s) and 


(O(s). 
Cx {i 1-0 = Fh 
Theorem 7. WU ( = oiseIs) = 3. 


Proof: Let W,(s)= Oll-a)Js]-& 


[A 
Then y' (S$) = Oo l-e)slu-a) = (0 [t-c)s] 


a ee 
and (ol = (0) =A. . Therefore YW and W have equal 


first moments, Since yi (3) ip" -a)s) (1-0) | a 
$20 [$20 


Af 


W, and y have finite second morents, 


y, Yrs) = Dlliralsut-a FL Wh-alsl) - a 


1-=a f= 


= F(ll-a)yisital = @ 


{[- 0. 


(ws) = KEY, (S)] 


Therefore |) and |) satisfy the same functional relationship with k(s). 
Therefore by theorem 6, W/s)=W,(s}, 


Now we can determine the relationship between H(u) and G(u). 


Theorem 8, {u) = = ie al- 0. 
ai é 
Proof: From theoren7,  (s) = 1PLl-G)s]-a ~ 
I—- 
By definition of generating functions {Ss} = Pred H/{u) 
0 
a [ve -a. |= ClicQ/SJ]-@ 
rrr. Y détu ora -Q 


Cog, 5 co ult-a)s 
tence 0 dui = ata | fre dl) - a] 


=)]@2 





U-eaen dHly} = \ e neues d Giv)- 
u/ O re, 


Let &(U) be defined as follows: f(vVizs G for UO 


g(vizQ& Fer upo 
of 00 
Then @ = { d€(u) = 4" dele. 
O™ 


Therefore a is the Laplace transform of £ bw), 


Hence j= e*” (1-a) id HI u) = in edGla)-) > a aici) 


Jo 
a 4 ei [li-a a)H(.)] = is “as” IG ate (ui) 
Therefore (j-Q)k(u] = G Geen) = ¢lu) 


and MM (u) = i | Gl ita} a] "or UFO 
Clearly H(u) = O and H(ee) = 1 as it should, 
To put the relationship of H(u) and G(u) in terms of random 
variables, we prove the following theorem. 
Theorem 9, The random variable w is distributed as the product of 
two indeperdent random variables, w, and w', where P(wy = 0) = a, 


P(w, = — = lea and w! has the distribution H(u). 


| | sy 
Proof: yes | = 0. ¢(l-a! Ce I< a Hlu} 


QO. + Ua) fe Vd H (yis~aal 
Pre dGlu 


ip (s) 


J 


fis 


ee 





So, P( wli-a) ¢uU|wro) = Hu). 
For later use a theorem due to Harris [10] is given with some 


remarks, Let {= logs ay © Logul seg) ) If qy 0, let VY=o. 

Theorem 10, If y<oo, Re(s) £0 , #0, then W(s) = Si t Mas} 
where M(s) is ee for s # 0, MU s) = M(s); Mj(s) = 0 (9) Sie 
M(-s) = M(s). 

Remarks: (a) Under the conditions of the theorem M(s) is real and 

poaitive when s is regal and negative, 


coy if E(2")<00, the r*" derivative of W(s) satisfies 


tr) = Ora iS| ~ 0 


W (s) 








(c) If Y= co, Y(s) and as many derivatives as exist approach 
zero exponentially as |s|+o, 

We will now determine some of the properties of G(u); however, 
it will be convenient to prove theorems for H(u) and then interpret 
the results in terms of G(u). Let h(u) = H' (u). 

Theorem 11, H(u) is absolutely continuous. 
Proof: Cramer [3] breaks up an arbitrary cumulative distribution 
function F(x) into three pieces: wile 
F(x) = Q) FAX) +a Fa ix) + G3 Fe ly) 
where ©,+0,+43 =1 Aiaz Oo a eS 
F, (x) is absolutely continuous. 


F(x) is a step-function and is equal to the sum of the saltuses 
of F(x) at all discontinuities which are less than or 
equal to x, 

F(x), the "singular" component, is a continuous function having, 


almost everywhere, a derivative equal to zero. 
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Lt X ; 
5) P, [+] = ~ 9 @ is the sum of an absosutely convergent 


trigonometric series and is thus an almost periodic function which 
comes as close to ay as we please for arbitrarily large values of t., 
By theorem 10, Lim W(t) =O, Therefore, for H(u), a, = 0 and H(u) 


is continuous, 


mn = 
Let h,,(¥) Fe { _ Wlit) dt ope iene 
Let dy = ot it Uz Wit) 
y= a ety duzd( We) 





mM .yn _it 
yer Er, d wilt) 


m aa 
(a) hu) = > sHoT | ylume y lume] + me Whe lt 


rn 


Examine each term of (A) as m =»©oo 


First iy (um) _— = ) te + Mos} a by theorem 10. 


ie 


-7 \ -LMvV 
| [mp Mlum) ¢ O(a] | e 
in| 
As we pass to the limit, this term is a function of continuous 
functions; hence, it must be continuous, Similarly the second term 


“imu “7 ] Ln 
‘leads to the expression ly (-emlé = | il Mim) + Oty) |e 


which converges to a ae funeeaa 


The last term in (A) is { 4 aaa — abled di. 


Since i) Fin de | < ) Fx) dv, it follows that 
_ | 











ea” 6 “eo 
Si t) ous 
: LA atl | dt S is car dt 
1) f | 
cS a are at 


a2 le 








} 
4174 
Km O| im) by corollary b to theorem 10 
for some constant k, 


okies  ) 
im} 
Hence this term approaches zero as m+ , Therefore, the continuous 
functions h(a) converge uniformly in some interval, ug Us U, 


to a continuous function h(u), 
26ED5, ally 


Hlu,) ~ H(u,) = lm " @ ~£ Wit) dt 


M~7 oO die aoe 





f° ee Ae 


a 


= lim " ou { , * Wit} dt 


M300 m 


= = lief hy, since h(a) converges uniformly 
to h(u) the limit sap be taken inside the integral and 
Hiv) = =f? h (vu) du 

Since H(u) is continuous at zero and we have shown ~bsolute con- 
tinuity elsewhere, it is absolutely continuous on an interval in- 
cluding the point zero, 

In case El 2K eco andr<%+K-1, integration by parts of (A) and 
corollary b to theorem 10 shows that the first r derivatives of 
h(u) are continuous if u # 0. -The integral expression for h(u) 
in terms of W(it) shows that ¥>r+/ implies n6?) (uy) is continuous 
at zero, This implies that G(u) =a +f gividv for u>0; g(v) 
is continuous forv #0, If E(2\)<0o, tnen e(T) (u) is continuous for u#0 


provided r¢y4+-} and is continuous for u = 0 provided rg 7-! . 


apis 








Corollary to theorem 1l, p_. = P(z, = r)40O uniformly in r, ril 


as N=>00, 


. -_ . ‘| = fe 
Proof: a - Gol) Gl 


+ [Gin] ~ 6, (Sr) 
Since G_(u) »G(u) uniformly for u>0O, G(u) is uniformly continuous 
for u>O and G(u) is right continuous at u = 0, the cegollary is 
proven, 
It is interesting to note that Yaglom [22] treated this subject 


in a manner similar to Harris, Yaglom defined 


Go, ea/- 
n=l Pho 


pf 


on 


Hlvi = Ply, €yly so} 


f 
Yn 


Hiy) = hm Hy ly) 


N-po 


oe) 
it 
lee] = je "d H ly) and proved the following theorem, 
Theorem 15, If mol , the variance of 24 is not zero, and £ (Zee 


then F LWit} (l-a}) ta) = (I-a) (me) ta 


st In this paper we followed the general line of Harris‘ argument 


Bs defined k(s) such that iy ( MMS) = *Mirabuls| ta] =a ; 


This led us to the conclusion that H(u) was the distribution function 


v4 
for (l-a)w which is the limit of (Pro) ar ; 
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MEAN NOT GREAT&R THAN ONE 

We have examined the case for the mean greater than one and 
now we turn to the situation which corresponds to the problem of 
subcritical or critical atomic piles, From theorem 3, the prob- 
ability of extinction is one and our fissioning process will 
eventually die out, 

Wher the mean is exactly one, the limiting distribution of the 
conditional random variable, y,, (which was defined in the last 
section to be the number of members of the nth generation, given 


that this number is not zero) is the same for any discrete branching 


process, 


Pnr 


Let a, = t=D ho ; Dee = a = ta = r | a #0) 


oo 
Kes) = a ce be the generating function of the p*. 


if the limit exists 


K = lim | J Flo 


a 
HAly) = PUY: <YIY, >0) 


Theorem 16, Cy eam c's) FO, c''[1)Zc0 then 


hm Hyly) = Hly) = /-e y =O 


Noo 
=O Y<O 


(The proof may be found in Yaglom [22].) 


Now we will consider the limiting distribution of F(s) in 


the case .2 ies 


ere 





Theorem 17, If K<! and C" £00 then 


him, Mar r = p* =| 
lim F (s) = | F(s) 
jroo % . 
ana FCF ISN= se Fis) + (leu) isl <| 
F{i}=| rij=eKk 


Eroet : 


= RIBRE poss ZatO 


ao b 
y=) ey ( Zn #0) . y definition of conditional 


probability 


— S S Elias) 


Since the summation is overr 0O 


2 |= mk 
— felsi= Foley} 
oi. F,10) 


The following theorem may be fourid in Konig fib}. If 8(s) is analytic 
in a domain, R, such that O(s,) = s, and | 9(s,) | = agi for some 


Sy in R, then, in this domain ESAS ease ee converges uniformly to 
O 


a continuous Pumietion D(s) such that HNL els] =a D (s) 


=| 


S.] =O Dis, 


SO 


25. 





f(s) is analytic in js|Z1. Since f(1) #1 and f'(1)<1, (Mé1), 


the s, of Konig's theorem is one, We can conclude that FA (Si f. 


. yi 


converges uniformly to a continuous function D(s) such that 


(2) D(f(s)) 2aD(s), D(1) =0, D(1) 21. 


» From 1), Fils) = | + | Els leary 





N=» N~—poo 


Lim F(s] = Els} = lim ‘ i (ata) 


(3) F(s) = 1 +-D(s) Kr by Konig's theorem, since convergence of 


F(s)— | 
mw” 


—- is uniform, the limit of Anhalt) exists. 
Zs 


From (3), then, D(s) = LE is)<i} kK 
From (2), D(F(s)) wr D{s) 


/ 


i 


| 


pM KC F(s}-1 
Ee Cis\y =] 4 pee by (3) 


| F(s)~1] 
Se KLE 


i 


1+ mFIs) ~u 


And, finally, FOS(s!} = mw Fis) $l!) 
From (3), F (i) | + De 


i 


A since P/i)=9_ by(2) 





One of the more interesting questions in the subject of branching 
processes is the determination of the number of generations to 
extinction when peé ibs 


let N=mintn: Z,,, 20] 


= oO 
Ol(s)=S N= | e> =~ ae NS 
no 
b, = Fn b, = I-p 
a = NS 
a j 
, iS bpre 
yy=0 


Theorem 18. The probabilities b, satisfy the recursion formula, 
b pose | on & (1-b,)} 


n+! 


Proof: b — Hx ji= Sp O\~ 
= “n+ Pn 42,0 | =. Pri Pato S] 


Theoren 19, The moment generating function of N is given by 


Ols) =} + (e°-1) 6, (5) 


(= 
Proof: @Ois)= Slp -~p jes 
eae ae as an f <0 
ar Nel, a) 
oe) oO 
Eas = : nS 
ae Mie i a 
REO n=0 | 
—_e > | Sins] ans ns nts 
= Es € t Treyoe = nso © j 
= | = 3 a" 
=! + > (e~- ] ae Freie 





So ©(s) is determined once 6,(s) is known, If we let 2 = e°, 





9,(s) is a power series whose coefficients are successive iterates 
, WY. oa : : ao ete 2 
of f*(b) 21 - f(l=b) since Diet = f(b) = £7, (bg), where £*(0) = 0 
a c=) 
£* (0) =4<1, (see Fatou[8].) A function of this sort is meromor- 
phic with poles ats2-=n logf@ , n21,2,. . « So 


\/ 


al 
Vr { : 
8,(s) =~ apie =; Hee ate Kap, 
i A Aare = fA 309 
converges everywhere except at the poles, (see Lattes [29] ) 
2 
is ae x;s1 is determined by X(xu 8) = f° (see Gx 
Yo is found by X(yo) = by = 1 = py or by using inverse functions 


x-l(f¥*(s)) s “x-l(s) and determining x-1(b,). 

If our process is of the type in which the "parent" does not 
die, then the total number of individuals is of interest, If p61, 
the probability that this number is finite is equal to one by 


theorem ha, 


Ly, 
j Sts rome | . 
Let S= > Zp; then PiS<a)= | for fg! 
Azo 
co 


= mec a See oles 
Let a, =P(S=r), gis)* > as for lat G ' 
Theorem 20. If m4 1, then g(s) = s f(g(s)). 
as) 
Proof: Let §$,= => 2» 7,20 for some n with probability one 
bey 
when 4 1 by theorem ha, Therefore, On converges to some random 


variable S with probability one. The random variable, S,, represents 


the total number of individuals in the generations from the zerouh to 


1% 
t 


the nth generation, Let $02 itn where 9/7 om,..) Since ae 
2 : | 
Let a Sis 2 Sen , where 4-- 52 and 3,,2 Ply =}. 


a a 
yW= 

aN 

= 


The event Y, = n can be divided into n mutually exclusive events 


with probabilities: 
P(z) = i)P(n-i individuals in the following k-1 generations) 
Louis clyc cer: 


=2 ee 





So ©(s) is determined once 6,(s) is known, If we let 2 = e®, 
91 (s) is a power series whose coefficients are successive iterates 
of f*(b) = 1 - f(1-b) since b,,, = £%(b,) = £7.,(b,), where £*(0) = 0 
f*'(0) Z4<1, (see Fatou!8].) A function of this sort is meromor- 
phic with poles at s=2-nlogm, n2i1,2,. . . So 

ais) = WY + fale, : LY, 

converges everywhere except at the poles. (see Lattes [29] ) 

X(s) _— x;s1 is determined by X( «4 8) = f*(X(s)),X (L)s Grae ae 
Yo is found by X(¥o) = bb SF 1 & Po or by using inverse functions 
x-l(f¥(s)) = “xel(s) and determining x-1(b,). 

If our process is of the type in which the "parent" does not 
die, then the total number of individuals is of interest. if #41, 
the probability that this number is finite is equal to one by 
theorem la, 

Let os in; then P(S<o) = for ME | 

nso we) 

Let Gq, = PS = ‘an gis)= zars for is] <?, 
Theorem 20. If wé1, then g(s) ss f(2(5)). 


7,=.) for some n with probability one 


Aa 
ri 
pt 
; 


Emoot: Tet Spas =. 


ey 


when 444 1 by theorem ha, Therefore, 5, converges to some random 
4 


variable S with probability one. The random variable, S,, represents 
the total number of individuals in the generations from the zero’ to 
‘4: 


° . 4 \J ee 3 =o f 
the nth generation. Let Se Where Wee oe oe, SLCC. — ee 
fee 
=< { \ —— cS a a“ 5 
Let 6. ES: =. y where 4, = 2 and 3, = er fea}. 


W=o 


The event Yy, = n can be divided into n mutually exclusive events 
with probabilities: 
P(2, = i)P(n=-i individuals in the following k-1l generations) 

fOre l= ee. cere 
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Starting with one individual at time zero, the probability of pro- 
ducing n individuals in the next k generations is the coefficient 


of gh 1391 G(s). 


Therefore P(y, = i)P(n-i in the next k-1 generations) 


‘ 

oo 
wat 
= 


: Ton 7 < 
ae . i 
=o. p; (coefficient eles Vas (3) )) 


a 


a» 


mae ° 
=‘ p,(coefficient of sin (sG,_4(s))*) 


in 
Jj 


coeiiicieny ote s 10 (3 ps (sG,_3(s))}+) 


- 
1 


Eh 


he 


coefficient of s™ in (7> ps (sGy_7(s))4) 


U6 


coefficient of s® in f(sG,_1(s)) 


» 


Therefore, G,(s) = - My nS" = f(sG,_1(s)). Clearly, g,(s) = sGp)(s) 
ne 0 
since P(y, =n) = P(Sy = n+ 1) and g.(s) = sf(g,_4(s)). Taking the 


limit as n approaches infinity, we see that g(s) = sf(g(s)). 


Using theorem 20, the expected total number of members of a 


; : ; Cee eo. an ! | 
family is: /(1)2. 2 92 fo/sil) = #1) gis}! + oe Tgisllg (sii 
“3 : ali “an te wy '@ \ me oF nt) t 
‘ins "$21 : '§a] 
~ r] r) ‘ 
sof Ft af) 
nv] 
need ohn f 
Then, a ics eat 
fy. 
‘| \ 
VY C ma oe —— Cui & Al t~ La] e ¥ 
eeiariy, 2 = —~———-/— 7 — + and the variance of S is: 
a ~ fA 
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ESTIMATION 
In the introduction we mentioned that the parameters, Pry must 
be empiraca ly determined for the given amount and arrangement of 
atomic material, The method of maximum likelihood can be applied; 
however, it seems necessary to assume that certain random variables 
are observable. 
Let Zny be the number of individuals in the m*) generation 


who have exactly r descendents in the (m+ 1)8% generation, Then 


=) 22 a 
7 Zz =a — g 
lezen = > V2.0 let S, => 2, 

r=0 y=0 r= 


Theorem 21. The maximum likelihood estimates of p,; and Va based on 


observed values of z,, for mén (n+1 observed generations) are 


me=0 - 
Proof: In order to construct the likelihood function, first we 


determine the conditional distribution of z2,,, r=0,l, .. . given 


Zn, 0 From our initial assumptions, each member of the me peneration 


has the same probability distribution associated with his progeny. 


Therefore the probability of k, members having no descendents is 


Ke { Z 
P, ‘ | where f) is the binomial coefficient and is equal to 
Ss STO o 


Va 


(Zn Ke)!ko! » Of the remaining 2, - ky members, say Zp - kK) have 


e ° ° {Zmn7 X ss e e 
one descendent, This probability must be | Pr » Continuing 
\ NG J ° b 


in this manner we get the multinomial functions: 


P(2mos Zl» ° ° °9 “mis 2 © ° n=l ,0> SgeD yl 9 0 © ey Atel sue ae) 


oS 
Ze 
Ce ee 
ee ceO | | 
Oe [Zmr)! 





Therefore P(2n95 21> + + + | 2n-1,0> 2n-l,ji> °° -)P(2n1,09 4n-1,1 


meh 2n-2,0 > 2-2, mee 7 Te “Pi19 > 213 + -+ 202) 
Zm- 
Yee 

m=0 TT, (Zme): 


Due to the Markov properties of na branching process this is the 
joint distribution of the 2,,,m=0,1,...,re0,l,... 


The only factor of the likelihood function which depends on the De is 


Se Za 
oS tm 
r=0 


CO n 
, the logarithm of this expression is > es Zmr | 0G, (p.3) 
rso 60M 


To maximize this exoression subject to the restraint, 3 F = | 
=0 


the method of Lagrangian multipliers may be used. 


oD Nn 
Let => = ae log Pr -/ p.} 
yao mM=0 
Then aL = i 
ory > [me a }-A) 
setting ol — ,we get S 2 yp forsT)=26,bore. 
DP, — mr ry 


ed 
But a p= 03 summing the above expression over r gives: 
r=0 Nn 
SS Zmr — ) 
rs0 m=o 
| The double sum of the 2... is merely S, So ) = 5, and 


n 
ie Zr 


in = S, 


18 


oo 
Ivey sins! pe recall that JA= 2." Pr and that VY Zip = 2me4s 
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Using these facts, it follows that: 
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— 

— m=o0 F=0 
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ION FLUCTUATION 

The number of ion pairs produced by a fast primary particle 
passing through an absorbing medium is a branching stochastic process 
(Moyal [18 | ). In this process the primary particle, which may be 
a gamma ray, has a certain charge, mass, and velocity or a given 
frequency; the thickness and atomic properties of the absorbing 
medium are known, The primary interacts with the atoms in the 
absorber, freeing electrons by the Compton effect, These electrons 
are the first generation and each may ionize another atom or atoms 
depending on their energy; secondary electrons may ionize more 
atoms, etc. The process ends when all the freed electrons are too 
Slow to produce further ionization, This results in some number of 
free electrons and, thereby, an equal number of ion pairs. 

In determining the distribution of this number we shall make 
the following assumptions; 

l. Successive ionizing collisions are statistically independent. 

2, The total energy loss of the primary is very much smaller 
than its initial energy and, hence, the loss of primary energy may 
be neglected, Therefore, we are not treating slow primaries or 
very thick absorbers, We also assume that energy losses due to 
radiation, nuclear interactions, etc. are negligible. 

3, The absorbing medium is homogeneous, 

4, The recording of data is delayed sufficiently to allow 
termination of the process. Photoelectric effects in the gas or 
chamber walls is neglected, 

Unlike the branching processes treated in the previous sections 
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the probability distribution of the number of progeny per individual 
in a given generation is dependent upon the generation, Thus one 
may expect different methods to be used, Moyal has applied the 
Laplace transform method to solve the Kolmogorov equations of a 
Markov process and this is outlined below, 
Let t be the thickness of the absorber 

N = the number of absorber atoms per unit volume 

6(E) = total cross section 

q = primary ionizing collision rate 

Q,, = P(n ion pairs are produced in a given collision) 

p_(t) = P( n ion ¢airs produced in a thiclgiges t) 


© 
M(s,t) = p,(tje"S , the generating function of p,(t) 
Nn=d 


Then (N 6(£) (dq) dt) is the probability of an energy loss between 
BK and —& + aE in an absorber of thickness dt. The primary ionization 


rate is Q= N So (E)dE 


oO 


The > G,= | and q(dt) = P( an ionizing collision in thickness Gai 
n=O : 
Since at least one ion pair is produced in a given collision, q, = 0. 


Based upon the above assumptions we have 
(4) Py(t, tta) = 2. py (ti) Pha) (ta) 


p, (St) = (1-9 (58) bn + (q15t))9,, + o(ft) where {,, is the 


Kronecker delta and 0(§t) goes to zero faster than (t as §t0. 


By assumption 1 and equation (A), 


M(s,titta) = M(s,t)M(s, a) 


Soe 
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OM(s,t} . by LS Mm Sibdie M\ (5,6) 
hs 


d §t 0 
_ len qe (M(s, 5-1) M(s,t) 
= fee, pelt alse) + 4560S 9,€™)-I Mis 
: Y=} 


a ws. q M(s,t) [> qe = 
o Ms) 
ot 


} oa M (s, ¢) [9 eS q (e"*~1} 


Comparing this equation with ay. = cy » it is clear that 


Mis,t) = exp[at > (e114 ] - o @RIs) 
n=) 


=e) 
where Q= as Ris)= > (e7"* 1) 4, 


n=! 


The inversion of the Laplace transform, M(s,t), is given by: 


C+ie 
Riv) 
Pr (Q)= si \ °° en du (see Widder [21] ) 
Core 


Evaluating this integral by the saddlepoint method (see Korn [15]), 


using the first term only, gives: 
eve 
PalQ) = [ar QR"s}]  &*P[Q(RISn)~ Sy R'(Sn)] 


where Sn is related to n by 
a2 KS 
! se nN 
n= -QR'(S,) =Q> (ke "4, 
K=\ 


PG. 





CO : 
Be tS,) = > Reape 


and c is a normalizing constant, 
(==) 
Since > PA(Q) =| 
Y= 


\ 


nO Me | l 
C= 2 |anQ R'(S,)] . eX Pp [Q(RISn) > Sn R (Sn) ] 


one 
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